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ABSTRACT 

We consider the classical self-dual Yang-Mills equation in 3-|-l-dimensional Minkowski 
space. We have found an exact solution, which describes scattering of n plane waves. In 
order to write the solution in a compact form, it is convenient to introduce a scattering 
operator T. It acts in the direct product of three linear spaces: 1) universal enveloping 
of su{N) Lie algebra, 2) n-dimensional vector space and 3) space of functions defined 
on the unit interval. 
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1 Introduction 

We consider the classical Yang-Mills field with value in su{N) algebra in 3+1-dimensional Minkowski 
space. We study the self-dual equation: 

Study of the self-dual Yang- Mills equation is important for understanding of QCD 
We fix the gauge following ^ : 

Ao-z = 0, A^+iy = 0, Ao+z = V2d^+iy^, A^_iy = V2do-z^. (1.2) 

Here Aq^z = Ao±Az, A^±iy = A^ =f iAy. 

The self-dual Yang- Mills equation for scalar su{N) field $ looks as follows 

□ $-2(7[9,+,,$,9o-.$] = 0. (1.3) 

It is associated with a cubic action flOl . Following P] we start looking for the solution of eq. (|1.3|) 



using perturbation theory in coupling constant g 

oo 

$(a;) = J2 (1.4) 

m=l 

Here ^^"^^ depends on the coupling as g'^~^. The first term satisfies a linear equation 

□ = 0. (1.5) 

We choose ^^^^ as a sum of n plane waves 

n 

$«(a;) = -z^T"^e-''=^"/(fc,). (1.6) 

Here T" are su{N) generators 

trT^T^ = S^K (1.7) 
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All kj are n different light-cone vectors k'j = and f{k) is a function with support on the light-cone. 
We will also use the following notations: 

x+iy {kj)o-z 

We have found explicit expression for The first two terms coincide with all other ^^'"^^ 

(m > 3) are different. 

Let us explain our solution. 

We shall use an abbreviation: 

<P{j) = T^^e-'^^^f{k,). (1.9) 
We introduce the following function: 

^<'"=Sw''" = /2^— = ^»'^^- <^-^»' 

Here /q is the modified Bessel function of the first kind. The integration contour is a circle around 
zero. We integrate in positive direction. 
Let us define a linear operator: 

^(ai,«2; Ji,j2) 

= g(i>{ji)P{juj2) ds e-'v(sa,g<P{ji)P{n,j2))v(sa2g<P{j2)P{ji,j2)). (l-ll) 

Here ji and j2 run through n values. Integration variables ai^2 belong to a unit interval [0, 1]. We 
shall consider T as an operator in direct product of n-dimensional vector space and as an integral 
operator. The kernel T(ai, 0^2; Ji, J2) takes values in an universal enveloping algebra of su{N). We 
are using -P(ji, ^2) which is defined by 

{Qh - Qj2y^ for Ji ^ h 
Pi3u32) = \ (1.12) 
for ji = ja- 

V 

The kernel T{a, a']j,j') depends only on the j-th and j'-th plane waves. It vanishes if j = j'. 
The function ( ITTlD is a kernel of an operator T 



iT){auji),(a2;j2) = T{ai,a2]ji,j2), 



whose index {a;j) takes a value in [0, 1] x {1,2, . . . ,n}. It acts on a "vector" {f)(a;j) which takes 
a value in the universal enveloping algebra as follows: 

iTf){a;j) = J2 daT{a,a';j,j'){f)(a';j'). (1.13) 

We call T the 'scattering operator'. 

We introduce two special "vector" s (see (pTOl)) 



(0)(«;i) = 0(j), (0o)(a;i) = 1- (1.14) 

— * — * 

For example, a scalar product of (f)Q and an arbitrary vector function / is equal to: 

Now all the notations are prepared in order to write down the solution of the self-dual equation 



T73|) which we have found: 

$fx) = -i$n ■ I ■ 

,1-r 

This is the main result of our paper. 



$(x) = -t$o ■ ( ) 0- (1-15) 



2 Solving recursion relations 

Substituting eq.( |1.4|) into the self-dual equation (|1.3| ) gives the following recursion relations for 

$('")(x): 

m—l 

□ $('")(x) =tgY: (9.+.,$(^'^9o^,$(™-^) - 9o_.$(^)9.+.,$ ("^-^■)). (2.1) 
It is convenient to introduce the following object: 

x+iy x+iy 

(2.2) 

It has the property: 

2A;i -^2 = X(l,2)(Qi -Qa). (2.3) 
The lowest component of $(x) in (|1.4| ) satisfies the free equation: 

= 0. (2.4) 



A solution is given by 

n 

= -z^T'^^e-^^^Vl'^i), ^] = ^- (2-5) 

We shall use the notation: 

4>{j) = T^^e-^''^^f{k^). (2.6) 



Let us analyze The recursion relation (2J.) is 



= ^9J2J2 {dx+iy(j){ji)do-z(l){32) - do-z(j){ji)d^+iy(j){j2)) 



n n 

= ^9Y. 0(ji)0(j2)X(ji,j2). (2.7) 

Here we have used the property X(ji, ji) = 0. 
Equation ( |2.7| ) can be solved by: 

In general, we can add any terms which satisfy the free equation to this solution $*^^)(x). We 
choose the minimal solution ( |2.8| ). 

Next, we consider an equation for ^^^\x). 

□<i.(3)(x) = z(7(9,.+,,<i>«9o-.$(2)-9o_.<i'(^)a.+,,$(2)) 

= ^^?'EE E 0(ji)0(j2)0(j3)(Q,.-g,3)"'(^(ji,j2) + x(ji,j3)) 

+ ^^?'E E Y.^i3im2m^){Qn-Qnr'{X{3i.3^)+X{32.3^)) 

= E 0(J1)0(J2)0(J3)(Q,,-Q,3)"'(^(J1,J2)+X(J1,J3)) 

+ w'T. E 0(ji)0(ji)0(j3)(Q,,-Q,3)"M^(ji,j3)) 



+ ^/E E <i>{ji)<PU2)4>{j2){Qn-Qnr'{X{ji,j2)) 

= ^^E E E Hji)Hh)Hh)iQn-Qn)-\Qn-Qjsr' 

+ ^/E E 0(ii)0(ii)0O3)(Q,i-g,3)-M^O-i,i3)) 
+ ^g'H E 0(ji)</'(j2)0(j2)(g,,-g,j-'(^(ji,j2)). 

Note that 

(%i + % + ^3)^ = ^(ii, J2)(gji - Qj^) + X{ji,j3){Qj^ - Qj^) + X{j2,jz){Qj^ - Qj,), 

{2kj, + kj,f = 2X{n,j,){Qj,-Qj,). 
Then the minimal solution of ^^^\x) is given by 

h j2(T^h) h{^j2) 
h j2{¥^h) 

-^/E E <t>{n)<t>{j2)(l>{32)\{Qn-Q32)-^- (2.9) 
Similarly, we have the minimal solution for ^^^\x): 

- 12 E E H 0(ji)0(j2)0(j3)0(j4)(g,,-g,j-'(g,.-g,3)"'(Q.3-Q.4)-' 

h 32 y^h) 33 {¥=32 ) 34, {¥"33 ) 

+ E E E '/'(ji)'/>(ji)'/>O2)0(j3)^(g,,-g,j-'(g,.-g,3)-' 

Jl 32{¥3l) 33{¥=32) 

+ E E E </'ai)</'O2)0(i2)0O3) 

32{¥h) 33{¥32) 

X {liQn - Qn)~\Q32 - Qn)-' + ^iQn - Q32)-\Qn - QnY^) 

+ E E E 0(ii)</'(j2)0(j3)0(j3)^(g,i-g,j-'(g,.-g,3)-' 

ii 32{¥3i) hi^h) 
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+ E E Hjimjimjimj2)hQn-Qn)-' 
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+ E E 0(ji)0(ji)0(j2)0(j2);.(g,.-Q, 



J2J 



+ E E 0(ji)0(j2)0(j2)0(j2)^(g,.-g, 



ill) 



(2.10) 



Let us now consider the general case We define coefficients D as follows: 

n n 

= E • • • E '/'O'l) • • • • • • , J-)- 

ii=i im=l 



(2.11) 



See fp.Gl). Then the recursion relation (12. Ih can be transformed into the recursion relations for D: 



(2.12) 



(/Cj, + ... + A;j„)2L>(ji,...,j„) 

m— 1 / i m 

E ^O'l' • • • >iO^(ii+l> • • • > J'm) E E ^Us,jt] 

1=1 \s=lt=l+l 



See ( |2.2|) . It is convenient to rearrange the sum ( |2.11|) into the following form: 



-^r-^EE E ••• E E---E 

x5.,+...+n.™0(jl)"^0(j2)"^ . . . 0(jz)"'/^(l"^2"^ . . . n). 



We introduce a notation: 

ni n2 ni 

D{r^2-- . . . n = D(jr'?^,sr?3, • • • , jT^), 

Xst = X {js , jt ) , Qst = Qjs- Qjt ■ 

The initial condition for m = 1 and the minimality of the solution are summarized by 



Then the recursion relations ( p.l2D are given by 
Y,nsntQstXstD{m''^ ...V 



s<t 



E pntXstD{r^2''^ . . . (s - l)"^-isP)D(s"^"P(s + l)"=+i . 

s<t p=l 



..I 



(2.13) 



(2.14) 
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t-1 nq 

+ 5^n,ntX,tD(^^2"^..(g-l)"''-lgf)D(g"«-P(g + l)"«+^..r) 

s<t q=s+l p=l 
nt — 1 

+ E E - p)XstD{r'2''^ ...{t- lY'-H'P)D{e'-'P{t + . . . T'). (2.15) 

s<t p=l 

The solutions D do not contain the factor Xst- The above recursion relation is decomposed into 
the set of recursion relations by comparing the coefficients at Xst {s < t). We have recursion 
relations: 

t-l riq 
q=s+l p=l 

Not all recursion relations ( |2.16| ) are independent. Note that 

t 

Q st ^ ^ Qu—i,u- 

u=s+l 



(2.16) 



All Qst can be expressed by Qu~i,u {u = 2, . . . ,1). The set of recursion relations ( p.l6|) are generated 
from 

Qi_i,tD(ri2"^..r) 

nt-i 

= y -^D{r'T"' . . . (t - 2)"'-2(t - if)D{{t - . . . r) 

^1 



(2.17) 



for t = 2, . . . ,1. But not all recursion relations ( p. 17 ) are needed to determine the coefficient D. 
With the initial condition ( |2.14| ), the recursion relation ( |2.17| ) for t = 2 is sufficient to obtain the 
answer. We must check that the solution satisfies the other recursion relation ( |2.17] ) for t = 3, . . . ,1. 
Combining eq.( |2.14| ) and eq.( |2.17| ) for t = 2, we have the recursion relation: 

D(ri2"2...ro (2.18) 
1 ' 



n2 — 1 



nin2 



Q^^ n2/^(r^~^2"2 . . . r ) + E ^i(^2 - p)/^(r^2P)D(2"2^P3"^ ...i 



no 



p=l 



To solve this relation, we first determine the coefficient D{1"'^2p). The simplest case is 

Dil2) = Q^,'Dil)Di2) = Q^l 

For > 1, we have 



D(12"^) = -Qr2 ^(12"^-') = ^(^12')"^- (2.19) 



Next we consider D{V^^2 



n2 ^ 



D(ri2"2) = ^— gi2^ (n2£'(i"^"^2"2) + niD(r^2"2-i)) (2.20) 



nir72 

Let us define C{ni,n2) by: 

Z}(l"i2"2) = ^^(Q-i)-i+-2-ic(^^, ,^2). (2.21) 
ni\n2\ 

Then the coefficient C{ni,n2) satisfies 

C(ni, na) = C(ni - 1, n^) + C(ni, n2 - 1). (2.22) 

With the initial conditions C{l,n2) = C{ni, 1) = 1, the solution is given by 

{n,+n2-2)\ 
^^"^'"^^=(n.-l)!(n2-l)!- 

Therefore Z)(l"i2"2) is determined as 

^nni2n2) _ (ni + n2 - 2)! ^ ^ 



Let us insert this relation into eq.(^TTB|). Then we obtain 



L)(l"i2"^..r) = —Q7^D(l'''-^2''\..r) 

rii 



+ E , /^ C(ni,p)(Qr2')"^''"-P(2".../"0- (2.24) 
ni! p! n2 



Using this relation recursively, we have 

D(ri2"2...r) = -^(Q72^rw(2''\..r) 

rii'. 



+ E %T^(E^(^i + l-^'P))(Qr2r^^'^(2".../"0. (2.25) 
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With the help of the relation, 



ni 



s=l 



C(ni,p+l) = ^C(.,p) 
we obtain the following recursion relations: 



(2.26) 



P2=0 



Til! P2! ^2 



ni\ (^2-^2)! ^2 

Here we replaced p2 with 77-2 — P2- 

Then, solutions of the recursion relations for D are finally given by 



(2.27) 



n2 nz ni_i l l 

P2 = lP3 = l p;_i = l S=l 



(pt_i + nt -pt - 1)! 



1 ^^pt_i+nt-pt 



(2.28) 



t=2 



with the convention pi = ni and pi = 1. We can see that this solution ( [^.28 ) indeed satisfies 
recursion relations ( |2.17| ). Combination of formulas ( |1.4| ), ( p.l3| ) and ( p.28| ) gives the perturbative 



solution of self-dual equation (|1.3| ). 

Some explicit examples of D are given by: 



D{l2...m) = Q-,iQ:^i...Q-^\^, 

Til- 



(2.29) 
(2.30) 



3 Scattering Operator 

Let us simplify the expression for the solution: 

00 

^(x) = E 

m=l 

00 m 00 CO 

= -^E^?"-'EE E ••• E E---E 

m=l 1=1 ji jiiy^ji) j,(^j;_i)ni=l ni=l 

00 

= -^E E E ^"^+-+"'-V(jir-..0ur^(i"^..ro. (3.1) 

'=1 {jiT^ji+l} {ni>l} 
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See ( ^.13 ) and ( 2.28| ). To obtain the last expression, we first exchange the order of summation 



with respect to m and 1. Then we use 5ni+...+n;,m in order to perform summation with respect to 
m. 

The coefficients D depend on the momenta by means of Q^^n = {Qji_i — Qji)^^'- 

D{rK..n = DirK..r;Q^l...,Qr\i). 

The factor j diverges for = jj. To regularize this divergence, we replace j in D with 
Pi-i^i = P{ji-i,ji). See eq.( p..l2| ). Then the solution of the self-dual equation can be written as 

oo 

$(x) = ^<l«(x). (3.2) 

1=1 

This is not the expansion in the coupling constant. is not homogeneous in g. Here 

n n n 

^^'nx) = E E---E^^'nji,---,jO, (3.3) 

ji=ii2=i ji=i 

and 

oo oo 
?ii=0 ni=0 

x0(jO"^+^ . . . ^{jir+'D{r^+' . . . r+'- Pu, Pi-i,i). (3.4) 

Especially for / = 1, 

|.(i)(ji) = -^^{J,) = -^T^ne-^^n-f{^k,). (3.5) 

We consider $'^^)(ji, j2)- 

^^'\3l.32) = -^9-' E E (#(J1))"^^(#(J2))"-+^ , ^ /7 1,^^^: . , {P,.r^-^^\ (3.6) 
nt^ont^o (ni + l)!ni!(n2 + l)!n2! 

Note that 0(ji) takes a value in the Lie algebra of the gauge group. The order of 0's needs to be 
considered. 

We shall use the integral representation for the Gamma function 

roo 

m\= ds e-' s"", (3.7) 
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to represent (ni + 722)! in (^.61). We shall obtain 



Now the summation in rti and ?t,2 factorizes. 

The sum in ni or n2 in the above equation can be expressed in terms of: 

Here V{a) is given by eg. ( P-.10| ). Note that Vi{a) has several representations: 

/df rl 
i/t+at ^ i^^2^y^ = Via). (3.10) 
/TTZ da 

Here Ji is the modified Bessel function of the first kind. The integration contour is a circle in the 
complex plane around zero in positive direction. 
Now the solution J2) can be expressed as 

^^'^\ji,j2) = -if dai ( cia2T(ai,a2; ji,j2)0(j2), (3.11) 
Jo Jo 

where T(ai, 0^2; ji, j2) is defined by 
T{ai,a2;ji,j2) 

= g<P{ji)P{juj2)l^ds e-'V(sa,g<P{J,)P{J^,J2))v(sa2g<P{J2)P{Jl,J2)). (3.12) 
Then <l>(^)(x) is given by 

$*^^)(x) = -z V V / dai da2T{ai,a2;jij2)(j){j2)- (3.13) 

The integration variable a was introduced by formula ( |3.9| ). 

Using a similar technique, we can write the coefficient D ( |2.28| ) as follows: 

n2 rii_i I 

E • • • E n X ^'"^-1.* 



P2=0 p,_i=0t=2 



(^2^12)"^ 'r}{stPt-i,tr-'''{st-,iPt,t+iy' {siPi 



X 



(ni + l)!ni! (nj + l)(pj(nt - pt)!)2 {m + 1)W 
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We replaced pt with pt + l and represented {pt-i + —PtV- by St-integration. In order to calculate 
we need to substitute this expression for D into the sum ( |3.4| ). Then we have 



E • • • E E • • • E n / rf^* e--^0(jOPi2 

„ . n ^ n « n « n + o 



ni=0 n;=0p2=0 p;_i=0t=2 



(^1 + \)\n\\ 



M K + (n, + l)!n,! ^''''^ 

Using formula (|3.9| ), we introduce integration over ai and ai and perform the summation in ni 
and n;: 

OO OO 712 i „l „l 

= E • • • E E • • • E n / e"''* / dai / rfa«5'0(ji)Pi2V'(aiS2fi'0(ji)Pi2) 

n2=0 n,_i=0p2=0 p;_i=0t=2"^0 JO 



X n ^^^^^'^^ '"^'I , ,v'vJ' : ^,,2 ' ^ ' ^ ViaiSig<l>im^,Mji)- (3.16) 



In order to calculate summation in rit and for j = 2, — 1, we need to calculate the sum: 

A (gf#(j,)P,-i,,)"'"P^#(j,)nf+i(gf+l#(jt)^M+l)^^ .o^7^ 

We first exchange the order of summation then replace rit with nt + pt'- 

pt=Ont=pt 



^ ^ {nt + l)ipt\{nt-pt)\y 



OO OO 



^p ^p ist9(i)ijt)Pt-i,tT'9(l)i3t)Pt,t+iist+ig(l)ijt)Pt,t+iy' ion 



Let us represent (n^ + + 1) -"^ by 



^ (3.19) 



nt + + 1 Jo 

This formula introduces aj-integration. The sum in (|3.17|) can be written as 



i"- £ ^p-):: £ Kwp^) , ,3.0, 



12 



By means of ( p..lUD the sum ( ^.17] ) is equivalent to 

/ datV{atStg(p{jt)Pt-i,t)g4>{jt)Pt,t+iV{atSt+ig4){jt)Pt,t+i)- 
Jo 

Then, we have 

rl rl rl J' rOO 

= -i dai da2... dai TT / rfste"^* 5'0(ji)-Pi2^(aiS25'0(ji)-Pi2) 

Jo Jo Jo t-oJo 



t=2 

l-l 

X 

t=2 



]j(v{atStg(j){jt)Pt~i,t)9(j){jt)Pt,t+iV{atSt+ig(j){jt)Pt,t+i)^^ 

t=2 

i / dai / da2 ■ ■ ■ dai 
Jo Jo Jo 

l[[J dste-''g(p{]t-i)Pt-i,tV{at.iStg<p{]t-i)Pt-i,t)V{atStg<p{^^^ 



X 

t=2 



Recalhng the definition of T (|3.12|) , we have 

^^^\3id2, ■ ■ ■ di) = -if dai f da2... [ dai 

Jo Jo Jo 

xT(ai,a2; ji,j2)T(a2,a3; j2,j3) • ■ - Tiai^i, af, ji^i, ji)(f){ji). (3.21) 

Therefore, we have 

^^^\x) = -«y"y" dai da2... dai 

xT(ai,a2; ji,j2)T(a2,a3; J2,J3) • ■ -Tiai^i, af, ji^i, ji)(f){ji). (3.22) 
See ( |3.3|) . Using the 'scattering operator', we obtain a simple representation for <l>^'\ 

¥^\x) = -i$o-{fy-% (3.23) 
See ( |1.14|) . Thus, we proved our main formula ( |1.15| ): 



$(x) = -200 ■ ( 1 + Eity] = -¥o • ( ) 0- (3.24) 



=1 



1 -T, 



4 Summary 

Introduction of an auxiliary linear space is typical for the solution of two-dimensional completely 
integrable differential equations ITT], In this sense, our solution is represented similar to 

solutions of two-dimensional classical completely integrable differential equations. 
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